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Abstract: In this paper, we explore the cosmological imphcations of interacting dark energy 
model in a torsion based gravity namely f(T). Assuming dark energy interacts with dark mat- 
ter and radiation components, we examine the stability of this model by choosing different forms 
of interaction terms. We consider three different forms of dark energy: cosmological constant, 
quintessence and phantom energy. We then obtain several attractor solutions for each dark energy 
model interacting with other components. This model successfully explains the coincidence problem 
via the interacting dark energy scenario. 
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I. INTRODUCTION 

Its now well- accepted in astrophysics community that 
the observable universe is in a phase of rapid expansion 
whose rate of expansion is increasing, so called 'acceler- 
ated expansion'. This conclusion has been drawn by nu- 
merous recent cosmological and astrophysical data find- 
ings of supernovae SNe la cosmic microwave back- 
ground radiations via WMAP [2| galaxy redshift surveys 
via SDSS [3] and galactic X-ray [j]. This phenomenon is 
commonly termed 'dark energy' (DE) in the literature, 
and suggests that a cosmic dark fiuid possessing negative 
pressure and positive energy density. Although the phe- 
nomenon of dark energy in cosmic history is very recent 
z ~ 0.7, it has opened new areas in cosmology research. 
Two important problems like 'fine tuning' and 'cosmic 
coincidence' are related to dark energy. It is thought 
that the most elegant solution to DE paradigm is the 
Einstein's cosmological constant [5] but it cannot resolve 
the two problems mentioned above. Hence cosmologists 
looked for other theoretical models by considering the dy- 
namic nature of dark energy like quintessence scalar field 
[5] , a phantom energy field [7] and f-essence [8] . Another 
interesting set of proposals to DE puzzle is the 'modified 
gravity' which was proposed after the failure of general 
relativity (GR). This new set of gravity theories passes 
several solar system and astrophysical tests successfully 

0. 

For the past few years, models based on DE inter- 
acting with dark matter or any other exotic component 
have gained great impetus. Such interacting DE models 
can successfully explain numerous cosmological puzzles 
including dynamic DE, phantom crossing, cosmic coinci- 
dence and cosmic age [10] and also in good compatibility 
with the astrophysical observations of cosmic microwave 
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background, supernova type la, baryonic acoustic oscil- 
lations and galaxy redshift surveys []Jj]. There are some 
criticisms on interacting DE models for not being favored 
from observations and that the usual ACDM model is 
favorable However the thermal properties of this 

model in various gravities have been discussed in litera- 
ture [isl . The model in which dark energy interacts with 
two different fluids has been investigated in literature. 
In [14], the two fluids were dark matter and another was 
unspecifled. However, in another investigation [iS^, the 
third component was taken as radiation to address the 
cosmic-triple-coincidence problem and study the general- 
ized second law of thermodynamics. In a recent investi- 
gation, the authors investigated the coincidence problem 
in loop quantum gravity with triple interacting fluids in- 
cluding DE, dark matter and unparticle [16]. 

There is no need to construct a gravitational theory 
on a Riemannian manifold. A manifold can be divided 
into two separate but connected parts; one with Rieman- 
nian structure with a deflnite metric and another part, 
with a non-Riemannian structure and with torsion or 
non-metricity. That part which has the zero Rieman- 
nian tensor but has non zero torsion is based on a tetrad 
basis, and deflnes a Weitzenbock spacetime. /(T) gravity 
is an alternative theory for GR, deflned on the Weitzen- 
bock non-Riemannian manifold, working only with tor- 
sion. This model flrstly proposed by Einstein for unifying 
the electromagnetism and the gravity. If /(T) = T, this 
theory is called teleparallel gravity [13, |l8| . It has been 
shown that with linear T, this model has many common 
features like GR and is in good agreement with some 
standard tests of the GR in solar system [l^- But in- 
troducing a general model /(T) backs to few years [l9[. 
This model has many features, for example Birkhoff's 
theorem has been studied in this gravity [20]. Earlier 
the authors in [21] investigated perturbation in f(T) and 
found that the perturbation in /(T) gravity grows slower 
than that in Einstein general relativity. Bamba et al [22^ 
studied the evolution of equation of state parameter and 
,^iantom crossing in /(T) model. Emergent universes 
eHii chameleon /(T) model is investigated in [23]. As a 
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thermodynamical view, there are many strange features. 
It has been shown in /(T) gravity, the famous formula 
of entropy-area of the black hole thermodynamics is not 
valid [25]. The reason goes back to the violation of the 
local lorentz invariance of this theory L25i, i26j . It shows 
that it is not possible to use the Wald conjecture [l^ 
for calculating the entropy as a Noether charge in /(T). 
The Hamiltonian formulation of f{T) gravity has been 
studied in [2^ and shown that there are five degrees of 
freedom. In this paper we study the triple coincidence 
problem: why we happen to live during this special epoch 



when Pa ~ Pn 



[28]. We investigate this problem 



in the framework of /(T) gravity. 

We follow the plan: In section II we introduce the basic 
equations as an autonomous dynamical system. In sec- 
tion III we choose a model for f{T) gravity. In section 
IV we working on numerical analysis of the stability and 
the evolution of the functions of the model in details. We 
conclude and summarize in section V. 



II. BASIC EQUATIONS 

One suitable form of action for /(T) gravity in 
Weitzenbock spacetime is [19] 

5 = ^ y d^x^iT + f{T) + Lm) 

Here e = (iet(e^), = SttG and is the tetrad (vier- 
bein) basis. The dynamical quantity of the model is the 
scalar torsion T and Lm is the matter Lagrangian. We 
start with the Friedmann equation in this form of the 
/(T) model [l9| 



(1) 



where p = Pm + Pd + Pr, while p^, pd and pr represent the 
energy densities of matter, dark energy and the radiation 
respectively. 

Another FRW equation is 
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(2) 



For a spatially flat universe (fc = 0), the total energy 
conservation equation is 



(3) 



where H is the Hubble parameter, p is the total energy 
density and p is the total pressure of the background 
fiuid. 

We assume a three component fiuid containing matter, 
dark energy and radiation having an interaction. The 
corresponding continuity equations are [l6| 



pd^^H{pd^Pd) = Ti, 

Pm + ^Hpm = 
Pr^3H{pr^Pr) = Ts, 



(4) 



which satisfy collectively (j3j) such that Fi + r2 + Fa = 0. 
We define dimensionless density parameters via 
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The continuity equations (j4]) in dimensionless variables 
reduce to 



dN = H—1 



^y^Z^ WdX + WrZ 



(6) 



-3x(l + Wd) ■ 
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dy 
dN 



dz 
dN 



-\- y -\- Z -\- WdX + WrZ 



-3y 



fr + 2TfTT 
rr2, 
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+ y + Z + WdX + WrZ 

fr + 2TfTT 



-3z{l + Wr) + 



where = Ina, is called the e-folding parameter. The 
coupling functions F^, i = 1,2,3 are in general functions 
of the energy densities and the Hubble parameter i.e. 
Vi{Hpi). The system of equations in (j6j) is analyzed by 
first equating them to zero to obtain the critical points. 
Next we perturb equations up to first order about the 
critical points and check their stability. Below for com- 
putation, we shall assume Wm = Wr = \ and Wd to 
be a general non-zero but negative parameter. We are 
interested in stable critical points (i.e. those points for 
which all eigenvalues of Jacobian matrix are negative) as 
these are attractor solutions of the dynamical system. 



III. /(T) MODEL 

To avoid analytic and computation problems, we 
choose a suitable /(T) expression which contains a con- 
stant, linear and a non-linear form of torsion, specifically 



/(r) = 2CiA 



(7) 



where a, Ci and C2 are arbitrary constants. The first 
and the third terms (excluding the middle term) has cor- 
respondence with the cosmological constant EoS in /(T) 
gravity [29|. There are many kinds of such these mod- 
els, reconstructed from different kinds of the dark en- 
ergy models. For example this form (7) may be inspired 
from a model for dark energy from proposed form of the 
Veneziano ghost [30^. But the linear term is needed to 
show the differences between our results in f{T) gravity 
from the Einstein gravity. Here we choose this model to 
simplify our numerical computations and for easier dis- 
cussion on the difference of our results with the same 
results in GR. It is the minimum model, but our equa- 
tions have been written for a general f(T) action. It 
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is possible by repeating the numerical steps as we done 
in this paper, discuss the stability of other models. Re- 
cently Capozziello et al [31] inviestigated the cosmog- 
raphy of f{T) cosmology by using data of BAO, Su- 
pernovae la and WMAP. Following their interesting re- 
sults, we notice that if we choose C2 = 0, a = and 
Ci = \^Ho{ftrnO — 1)7 than we can estimate the parame- 
ters of our proposed f{T) model as a function of Hubble 
parameter Hq and the cosmographic parameters and the 
value of matter density parameter. 

IV. ANALYSIS OF STABILITY IN PHASE 
SPACE 

In this section, we will construct four models by choos- 
ing different coupling forms F^ and analyze the stability 
of the corresponding dynamical systems about the crit- 
ical points. We shall plot the phase and evolutionary 
diagrams accordingly. For this reason, we must find the 
critical points of the ([6]), and then we linearize the system 
near the critical points up to first order. 




FIG. 1: Model I: Phase space for Wd = -1.2, b = 0.5, a = 0.5. 
It shows an attractor behavior. 



A. Interacting model - I 

We consider the model with the following interaction 
terms 

ri = -6bHpd, T2 = Ts = 3bHpd, (8) 

where 6 is a coupling parameter and we assume it to be 
a positive real number of order unity. Thus (|8]) says that 
both matter and radiation have increase in energy den- 
sity with time while dark energy loses its energy density. 
Therefore it is a decay of dark energy into matter and 
radiation. 

Using ([8j), the system ([6]) takes the form 

—— = -3x(l + ^rf) + 3x — -66x, 

dN V 1 + a / 

^ = -3^ + 3, ( ) + 36x, (9) 

^ = -3.(1 + + 3.( " + ^ + W^'f + "-^ ) + 36.. 

The critical points for this model are obtained by equat- 
ing the left hand sides of (|9|) to zero. We obtain four 
critical points: 

• Point Ai: ((i±^)li±^^ 

• Point Bi : 

• Point Ci : 

• Point 1^1 : (0,0,|(l + a)(|-6)) 

The eigenvalues of the Jacobian matrix for these criti- 
cal points are: 



• Point : Al = 3(l+^d+26), A2 = 3(^^+36), A3 = 
-l + 3(^d + 36), 

• Point Bi : \i = 3{b - 1),A2 = -3(6 - 1),A3 = 
-3(l + 26 + ^rf), 

• Point Ci : Al = -1, A2 = 3(1-6), A3 = -3(^^+36), 

• Point Di : Al = 1, A2 = 4 - 36, A3 = 1 - 3wd - 96 

Point Al is stable when one of these conditions is satis- 
fied: 

Wd < -3, 6 < 1/18. (10) 
Wd > -3, Wd < -10/9, 6 < 1/18. (11) 

Wd>-Y^ m<0, 6<-i(l + ^rf). (12) 

Bi is an unstable critical point since if 6 > 1 than Ai > 
but A2 < 0. Ci is stable if 6 > l,!^;^^ > —36. Similarly 
Di is unstable since Ai > 0. 

In figures (1-6), we plot the parameters of model-I. In 
figure 1, we observe that the dimensionless density pa- 
rameters evolve from their currently observed values to 
vanishing densities. In figure- 2, we observe similar behav- 
ior when plotted against e- folding parameter: the density 
parameters evolve from their current values to zero. In 
figure- 1, we deal with phantom energy, in figure-3, the 
quintessence case while in figure-5, the cosmological con- 
stant case. 



We study another model with the choice of the inter- 
action terms 

Fi = -SbHpd, F2 = 3bH{pd-pm). = SbHpm. (13) 



(0,0,0), 

(0,(l-6)(l + a),0), ^' Interacting model - II 
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FIG. 2: Model I: variation of x, z as a function of the N 
ln(a). The initial conditions chosen are x(0) = 0.7, y(0) 
0.3, z{0) = 0.01, Wd = -1.2 and b = 0.5. 



FIG. 4: Model I: variation of x,y,z as a function of the N ■ 
ln(a). The initial conditions chosen are x(0) = 0.7, y(0) : 
0.3, z{0) = 0.01, Wd = -0.5 and b = 0.5. 
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FIG. 3: Model I: Phase space for Wd = -0.5, b = 0.5, a = 0.5. FIG. 5: Model I: Phase space for Wd = -1,6 = 0.5, a = 0.5. 



This model effectively describes the situation when dark 
energy loses energy density to matter while the radiation 
density increases due to interaction with the matter. 



dy 
dN 

dz 
dN 



X -\- y -\- Z -\- WdX + WrZ ^ 

1 + a 
?>hx — 3x(l + Wd)^ 

X -\- y -\- Z -\- WdX + WrZ^ 



1 + a 
+36(x -y)- 3y, 



= 3z(^ 



X -\- y -\- Z -\- WdX + WrZ ' 



1 + a 
-\-3by — 3z{l + Wr)^ 



(14) 



(15) 



(16) 



• Point A2 

• Point B2 

• Point C2 



(0,0,0), 
(0,0,1 + a), 

(0,(l-36)(l + a),36(l + a) 



• Point D2 : (- 3(2+^a)(6-^;d-l)(6-^;d-7/3)(a+l) 



3 K;d3 + (16-3 6)'w;d2 + (_6 6+27)^^ + 14+62+6 5 
b3(a+l)(6-Wd-7/3)(6-Wd-l) 



3 Wd3 + (16-3 6)t(;d2 + (_6 6+27)^^ + 14+62+6 

3(a+l)(6-Wd-l)6^ N 

' 14+16 Wdb-\-b^-3 Wd^b-\-3 Wd^-\-27 Wd+6 / 

The eigenvalues of the Jacobian matrix for these critical 
points are: 

• Point A2 : Xi = -4, A2 = -3(1 + 6), A3 = 3(1 + 

Wd - b), 



There are four critical points: 



Point ^2 : Ai = 4, A2 = 1 - 36, A3 = 7 + 3{wd - h), 
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FIG. 6: Model I: variation of x^y^z as a function of the N = 
ln(a). The initial conditions chosen are x(0) = 0.7, y(0) = 
0.3, z(0) 0.01, Wd -1 and b 0.5. 




FIG. 7: Model II: Phase space for wa — —1.2,6 = 0.5, a = 
0.5. 

• Point C2 : Ai = 3(1 + 6), A2 = 36-1, A3 = 3(^^ + 2), 

• Point L>2 : Ai = -3(^^ + 2), A2 = 3(6-^^-1), A3 = 
-7-3(6-^^) 

A2, D2 are conditionally stable if 6 > l^Wd (for A2) and 
Wd > —2 and b < 1 -\- Wd (for D2) . But B2 and C2 are 
unstable since Ai > 0. 

In figures (7-12), we show the dynamics of Model-II. 
Attractor solutions are shown in figures 7,9 and 11. In 
figure 8, we see that the energy density of dark energy 
decays like quintessence. Also the energy density of ra- 
diation first increases till N ^ 1.6 and then starts de- 
creasing. The matter density always decreases and ap- 
proaches zero nearly N ^ 3. In figure 10, the energy 
density of dark energy increases rapidly behaving like 
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N 

FIG. 8: Model II: variation of x, 2/, 2 as a function of the N = 
ln(a). The initial conditions chosen are x(0) = 0.7, y(0) = 
0.3, z{0) = 0.01, Wd = -1.2 and b = 0.5. 




FIG. 9: Model II: Phase space for Wd = -0.5, b = 0.5, a = 
0.5. 



phantom energy, energy density of matter rises almost 
exponentially at later times, while radiation density in- 
creases slower compared to both matter and dark energy. 
These novel behaviors appear on account of interaction 
between three components. In figure 12, the dark en- 
ergy density behaves like quintessence, while matter and 
radiation density falls with expansion. 

C. Interacting Model - III 

Let us take the interaction terms [l6| 

T^ = -6bi^H-^PdPr, T2 = rs = 3bK^H-^PdPr- (17) 



0.1 0.2 0.3 0.4 0.5 



2 4 6 8 10 
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FIG. 10: Model II: variation of x, y, z as a function of the N — 
ln(a). The initial conditions chosen are x(0) = 0.7, y(0) = 
0.3, z(0) = 0.01, Wd = -0.5 and b = 0.5. 



FIG. 12: Model II: variation of x,y,z as a function of 
the N = In (a). The initial conditions chosen are x(0) = 
0.7, y{0) = 0.3, z{0) = 0.01, Wd = -1 and b = 0.5. 




FIG. 11: Model II: Phase space for Wd = -1,6 = 0.5, a = 0.5. FIG. 13: Model III: Phase space for Wd = -1.2, b = 0.5, a 

0.5. 



The system in ([6]) takes the form 
dx 



Point Ds : (1 + a, 0,0), 



( 



dN 
dz 
dN 



3z 



X - 






- WdX - 


- 






1- 


f a 




xH 






-WdX^ 








1 - 


fa 








-z^ 


- WdX H 





• Point E'^ ' _2il±EA _l±wd) 

^^^^^^ ^3-^96' 96 ' 66 /' 



^ — 3x — 3wdX — ISbxz^ 

)-3y^9bxz, (18) •Point Fa: ( ^^61+67^2 



1 + a 

There are six critical points: 

• Point As : (0,0,0), 

• Point Bs : (0,1 + a, 0), 

• Point Cs : (0,0,1 + a), 



(wd+6 6Q;+6 6-|) 
6 6- 

3 6q! -3 6+36 6^q;+18 6^Q!^+9 Wd6Q!+18 6^+Wd^-2 ^+9 ^^6+^ 



■ 3z — 3wrZ + 96xz. 



36(6 6Q;+6 6-i+2 Wd) 



Wd(36a+36-^+Wd) 
36(6 6Q!+6 6-f +21 



The eigenvalues of the Jacobian matrix for these criti- 
cal points are: 

• Point As: Xi = -3, A2 = -4, A3 = -3(1 + Wd), 

• Point ^3 : Ai = 3, A2 -1, A3 = -3wd, 
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FIG. 14: Model III: variation of x,y,z as a function of 
the N = In (a). The initial conditions chosen are x(0) = 
0.7, y{0) = 0.3, z{0) = 0.01, Wd = -1.2 and b = 0.5. 



FIG. 16: Model III: variation of x,y,z as a function of 
the N = In (a). The initial conditions chosen are x(0) = 
0.7, y{0) = 0.3, z{0) = 0.01, Wd = -0.5 and b = 0.5. 
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FIG. 15: Model III: Phase space for Wd 
0.5. 



-0.5,6 = 0.5, a = FIG. 17: Model III: Phase space for Wd = -1,6 = 0.5, a 
0.5. 



• Point C3 : Ai = Sw^^ X2 = 3(1 + w^)^ A3 = 9 6a + 
96-l + 3^rf, 

• Point : Ai = 4, A2 = 1 - 96 - 96a, A3 -96a - 
96 + 1 - 3wd. 

As is stable for Wd > —1. ^3 , 1^3 are unstable. Cs is 
conditionally stable when Wd < —l^ b < ^pq^. 

In figures (13-18), we have plotted the cosmological pa- 
rameters of model-Ill. In figures 13, 15 and 17, we show 
the attractor solutions of the differential equations equa- 
tion. In figure 14, we observe the oscillatory behavior 
of dark energy and other cosmic components. It shows 
that when DE energy density decays then correspond- 
ing densities of dark matter and radiation increases and 



vice versa. Figure 16 shows that all forms of energy den- 
sities vanish by N ~ 4.5 From figure 18, the radiation 
density stays zero while matter energy density decreases 
and vanish by ~ 2. The dark energy density increases 
by 2 while it decreases and stays constant at later 
epochs. 



D. Interacting Model - IV 



Consider another model with the interaction terms Ui 



Fi = -Sbhi'^H ^pdpr, 
F2 = 3bi^^H-\pdPr- 

Ts = Sbn'^H'^PmPr- 



PmPr 



(19) 
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The system in (|6]) takes the form 

3x{ ] — 6x — 3wdX — 9bxz, 

V 1 + a / 



FIG. 19: Model IV: Phase space for = -1.2, 6 = 0.5, a = 
0.5. 



dx 
dN 

dy ^ ( 

dN = H 



X - 






- WdX - 


- WrZ 






1- 


f a 




X H 






- WdX H 


- WrZ 






1 - 


fa 








-z^ 


- WdX H 


- WrZ 



) - 3^ + %{xz - 
^ — 3z — ?>WrZ + 96?/z. 



1 + a 

There are seven critical points: 

• Point ^4 : (0,0,0), 

• Point ^4 : (1 + a, 0,0), 

• Point Ca : (0,l + a,0), 

• Point : (0,0, 1 + a), 

• Point : (^\ 

• Point F4 : 

36(36+36a-f ) / ' 

• Point G4: (0,^,-^). 

For points A4, ^4, C4, 1^4, G4, the eigenvalues of the Ja- 




± 1 i+wd A 

3 6(1+K;d) '96' 3 6 ^ ' 



Wd+3b+3 6a-l/3 _ (- f +Wd) (wd+3 b+3 ba- f ) 
3b ' 3b(3b+3ba-f) 



FIG. 20: Model IV: variation of x,y,z as a function of 
the N = In (a). The initial conditions chosen are x(0) = 
0.7, y{0) = 0.3, z{0) = 0.01, Wd = -1.2 and 6 = 0.5. 



• Point G4 : Ai = —3wd^ 



cobian matrix 


are 


• Point A4 


: Ai 


• Point B4 


: Ai 


• Point C4 


: Ai 


• Point D4 


: Ai 



-3,A2 



-4, A3 



-3(1 + ^d). 



-3wd, A3 96(1 + a) - 1, 



-3wd-9b-9ba^l, 



. _ V3y'b(l+a)(3 ba+3 b/3a-4/9+3 b+ b) 
^2 — b(l+a) ' 

. _ VS^/bil-\-a){3 ba-\-3b/3a-4./9-\-3 b+ b) 
— b(l+a) 

It is observed that in this model, A4 is stable for 
Wd > —I ^ is stable for k;^^ < —1. C4 is unstable. 
D4 is stable for Wd < —1,6 > ' -"-^'^ possible 

to determine the stability of the point £^4. But G4 is 
unstable. 

In figures (19-23), we give description about model- 
IV. Figures 19,21,23 represent the phase space diagrams 
for different forms of dark energy. Figures 22 and 23 
shows that energy density of dark energy decreases while 
in figure 20, the density of DE increases. It's the typical 
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FIG. 21: Model IV: Phase space for Wd = -0.5, b = 0.5, a = 
0.5. 



In this paper we discussed the stabihty of the interac- 
tive models of the dark energy, matter and radiation in 
a FRW model, for a general f{T) theory. We derived 
the equations and show that why we have some attractor 
solutions for some specific forms of the interactions. We 





FIG. 22: Model IV: variation of x, y,z as a, function of the N ■ 
ln(a). The initial conditions chosen are x(0) = 0.7,2/(0) 
0.3, z{0) = 0.01, Wd = -0.5 and b = 0.5. 



FIG. 23: Model IV: Phase space for Wd = -1,6 = 0.5, a 
0.5. 



numerically integrate the equations and show that the 
evolution of the dark energy density mimics three diffract 
behaviors phantom, quintessence and cosmological con- 
stant in some interactive forms. We like to comment that 
this interaction is purely phenomenological required to 
meet some observational consequences. Since the phase 
space of the system of the evolutionary equations has 
a definitive end point, trackers are exist and depending 
on the initial conditions, there are different kinds of the 
trackers. 



behavior of the phantom fields. Thus, our model pre- 
dicts the correct evolutionary schem for the DE density 
in regime of phantom. 

V. CONCLUSION 

/(T) gravity is a powerful and novel theory for ex- 
planation of the acceleration expansion of the universe. 
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